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We calculate the two-photon 15-25' spectrum of an atomic hydrogen Bose-Einstein condensate in 
the regime where the cold collision frequency shift dominates the lineshape. WKB and static phase 
approximations are made to find the intensities for transitions from the condensate to motional 
eigenstates for 25 atoms. The excited state wave functions are found using a mean field potential 
which includes the effects of collisions with condensate atoms. Results agree well with experimental 
data. This formalism can be used to find condensate spectra for a wide range of excitation schemes. 



I. INTRODUCTION 



In the recent experimental observation of Bose-Einstein condensation (BEC) in atomic hydrogen the cold 
collision frequency shift in the \S-2S photoexcitation spectrum |^] signalled the presence of a condensate. The shift 
arises because electronic energy levels are perturbed due to interactions, or collisions, with neighboring atoms. In 
' the cold collision regime, the temperature is low enough that the s-wave scattering length, a, is much less than the 
c/j . thermal de Broglie wavelength, = y/h'^ /2TrmkBT, and only s- waves are involved in the collisions ||]. 
O ' The cold collision frequency shift has also been studied in the hyperfine spectrum of hydrogen in cryogenic masers 
c/J [Ql, and cesium and rubidium in atomic fountains. Theoretical explanations of these results and other work 
on the hydrogen IS-2S spectrum P,p^ have focused on the magnitude of the shift, as opposed to a lineshape. In this 
article we present a calculation of the hydrogen BEC 15-25' spectrum. We also describe how the formalism can be 
used for other atomic systems and experimental conditions. 



m 

^ ' A. The Experiment 

o 

I The experiment is described in |l|,H, and we summarize the important aspects here. Hydrogen atoms in the 15, 
OO F = I, mp ^ 1 state are confined in a magnetic trap and evaporatively cooled. The hydrogen condensate is observed 
in the temperature range 30-70 fiK and the condensate fraction never exceeds a few percent. Nevertheless, the peak 
, density in the normal cloud is almost two orders of magnitude lower than in the condensate and in this study we will 
neglect the presence of the noncondensed gas. 



, The two-photon transition to the metastable 25, F = 1, mp — 1 state (r = 122 ms) is driven by a 243 nm laser 



beam which passes through the sample and is retroreflected. In this configuration, an atom can absorb one photon 



■ from each direction. This results in Doppler-free excitation for which there is no momentum transferred to the atom 



and no Doppler-broadening of the resonance. An atom can also absorb two co-propagating photons and receive a 
momentum kick. This is Doppler-sensitive excitation, and the spectrum in this case is recoil shifted and Doppler- 
broadened. The photo-excitation rate is monitored by counting 122 nm fluorescence photons from the excited state. 
. ^ ' ^ typical laser pulse of 500 fis, fewer than 1 in IC* of the atoms are promoted to the 25 state. 25 atoms experience 

Jv>( , the same trapping potential as 15 atoms because the magnetic moment is the same for both states, neglecting small 
; I ■ relativistic corrections. 

[ The natural linewidth of the 15-25 transition is 1.3 Hz, but the experimental width, at low density and temperature, 
is limited by the laser coherence time. The narrowest observed spectra, obtained when studying a noncondensed gas, 
have widths of a few kHz |ll]. For the condensate, the cold collision frequency shift is as much as one MHz and it 
dominates the lineshape. 



B. Mean Field Description of the Spectrum 



The frequency shift in maser and fountain experiments has traditionally been described using the quantum Boltz- 
mann equation 0,^,1). In this picture, the frequency shift is the net result of the small collisional phase shifts 
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arising from forward scattering events in the gas. A mean field description, however, is more convenient for studying 
an inhomogeneous Bose-Einstein condensate. We will derive this picture in detail, but we summarize the results 
here. Collisions add a mean field energy to the atom's potential energy. For a 2S atom excited out of a conden- 
sate the mean field term is (5£'2s(r) — 47r?i^ai5_2s"'iS'(r)/TO. For a IS* condensate atom the mean field term is 
5Eis{v) = 47r?i^ais_isrii5(r)/m. (The fraction of excited 25* atoms is small, so 25-25 interactions can be neglected.) 
The ground state s-wave triplet scattering length has been calculated accurately {ais-is — 0.0648 nm The 
15-25 scattering length, however, is less well known {ais-2S = — 1.4± 0.3 nm from experiment |^ and -2.3 nm from 
theory ||). 

We denote the sum of the magnetic trap potential, V{r), and the mean field energy, 5Ex{r), as the effective potential, 
V^^^{r) (Fig. Here x is either 15 or 25. For 15 condensate atoms, the effective potential in the condensate is 
flat. Because ais~2S < and the condensate density is large, 25 atoms experience a stiff attractive potential in the 
condensate which supports many bound 25 motional states. 




radius [/im] 

FIG. 1. Effective potentials for 15 atoms in the condensate and excited 25 atoms. Selected single particle wave functions 
are displayed at the height corresponding to their energy. The dashed lines are the magnetic trapping potential V{r), which 
is identical for 15 and 25 atoms. The thin solid lines are the effective potentials, which include the mean field interaction 
energy. The vertical light solid lines indicate allowed Doppler-free transitions from the condensate, which must preserve mirror 
symmetry. The potentials and condensate wave function are for a peak condensate density of 5 x 10^^ cm~^ and a magnetic 
trap oscillation frequency of 4 kHz, which are characteristic conditions for a hydrogen BEG and a strong confinement axis of the 
trap [1,2]. The scattering lengths used in the calculations are ais-is = 0.0648 nm and ais-2S = —1.4 nm, and the chemical 
potential is ii/ks ~ 2 ^K. The 25 levels form a near continuum of motional states in an anisotropic three dimensional trap. 

The 15-25 spectrum consists of transitions from the condensate to 25 motional eigenstates of the effective 25 
potential. For Doppler-free excitation, the final states are bound in the BEC well. Doppler-sensitive excitation 
populates states which lie about ?i^fcg/2mfcB = 643 fjK above the bottom of the 25 potential, where Hko is the 
momentum carried by two laser photons. The latter states extend over a region much greater than the condensate. 
Because the excited levels are so different for Doppler-free and Doppler-sensitive excitation, we must treat the two 
spectra independently. 

The rest of this article presents a derivation of the effective potentials and a quantum mechanical calculation of the 
BEC 15-25 spectrum. 



II. 1S-2S PHOTOEXCITATION SPECTRUM OF A HYDROGEN BOSE-EINSTEIN CONDENSATE 
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A. Hamiltonian 



We start with the many-body Hamiltonian for a system with N atoms 



= ^ + Hj"* + V{rj) ) + H^"' + ff™", (1) 



3 = i 



where p^, r^, and are the momentmii operator, position operator, and internal state Hamiltonian respectively 
for particle j . V{r) is the magnetic trapping potential, which is the same for 15 and 25 atoms. 

7?'°" is the atom-laser interaction. After making the rotating wave approximation, it can be written 

jjias ^ J2 [(|25)(15|)je"*4"''* + (|15)(25|)je'4"'^*] , (2) 

where v is the frequency of the laser field {2hh' w E2S — Eis = Eis-2S on resonance). The laser beam is uniform 
over the condensate, so we treat the excitation as a standing wave consisting of two counter-propagating plane waves. 
The effective 2-photon Rabi frequency for Doppler-free excitation jl4j , 



is uniform in space. Here, / is the laser intensity in each direction, M2s,\s = 11-78 |l^] is a unitless constant, c is the 
speed of light in vacuum, a is the fine structure constant, and i?oo is the Rydberg constant. For Doppler-sensitive 
excitation, 

f^i?s(r) = f^i3s(e'"«^+e-''=«^), (4) 

where I^ds — ^ufI'^- 

jjcoii describes the effects of two-body elastic collisions. In the cold collision regime, the interaction can be repre- 
sented by a shape independent pseudopotential fl^ corresponding to a phase shift per collision of ka, where fik is the 
momentum of each of the colliding particles in the center of mass frame, 



^'°" = ^E'^('-'-^^-)K5-i5(|15)(15|).(|l^)(15|),+ai5-25(|e^Si)(e^S^|),, 

+a25-25(|25)(25|).(|25)(25|),]. (5) 



The sum is over N{N — l)/2 distinct pairwise interaction terms. The 15-25 interaction projection operator is written 
in terms of 



3 _ |15).|25), + |25).|15), 
6 ^uii'j — ID; 



because the doubly spin polarized atoms collide on the e"^E+ potential during s-wave collisions p^ . As mentioned 
above, the 25-25 scattering term is negligible for the hydrogen experiment, but it is included here for completeness. 

Inelastic collisions, such as collisions in which the hyperfine level of one or both of the colliding partners changes, 
will contribute additional shifts which arc not included in this formalism, but these effects are expected to be small 
in the experiment ||^. 

B. System before Laser Excitation 

We make the approximation that the system is at T = 0, and all atoms are initially in the condensate. T — 
models have accurately described many condensates properties [ pT| , and we leave finite temperature effects for future 
study. The state vector can be written 

|vl/o) = |15^0_^,0). (7) 

N terms 
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where \1S, 0) refers to the single particle electronic and motional state of an atom in a IS condensate with N atoms. 
We use the ket notation (|a; 6; ..; c)), in which the entry in the first slot is the state of atom 1, the second entry is the 
state of atom 2, etc. 

Minimization of (5'o|i?|5'o) leads to the Gross-Pitaevskii, or nonlinear Schrodinger equation |l8yi9|] for the single 
particle BEG wave function, ■0(r) — (r|0), 

A^^(r) = + K/^(r)) V(r). (8) 

The effective potential is V-^g'^ {r) — V{r) + Un{r), where C/ — 47r?i^ai5_is/m. Here, n{r) — N\tp{r)\'^ is the density 
distribution in the N-particle condensate. One can interpret |?/;(ri)p as the probability of finding condensate particle 
i at position r^. 

The kinetic energy is small and can be neglected. This yields the Thomas-Fermi wave function ||20[| , 

^Pir)^{ in{0)~V{r)/u]'^^ V{r)<n{0)U ^ (g) 

otherwise 

where n{0) is the peak density. The density profile is the inverted image of the trapping potential. The chemical 
potential is /u(iV) = Un{0), and it is equal to V^g^ inside the condensate. The energy of the system before laser 
excitation is the minimum of {'$q\H\^q). It satisfies n{N) = dEo/dN and is given by 

Eo = ^Nfi{N). (10) 
From now on, when writing /i we will drop the explicit dependence on N. For a cylindrically symmetric harmonic 

/ 3 3/2 \'^/'^ 

trap, it can be shown that n(0) = I IbNui^w^Wz /fi o-is-is ) /^tt, where Wr and Wz are the angular frequencies for 
radial and axial oscillations in the trap. 




C. System after Laser Excitation 

To describe the system after laser excitation we must find the orthonormal basis of 25 motional wave functions and 
their energies. This is done by minimizing where 



1$,,,) =5|25,z;...;25,z;15,0;...;15,0) 



(11) 



q terms 



N— q terms 



is a state with q 2S atoms in 25* motional level i. The operator S symmetrizes with respect to particle label. We will 
show below that the state vector of the system after laser excitation is actually expressed as a superposition of such 
terms, but for now we need only consider a single 

Calculating {^q^i\H\^q^i) involves a somewhat lengthy calculation. Details are given in appendix ^ and the result 

is 



{^q,^\H\^q,^) = E'o 



-q (25, i 



E'„+q{Eis-2S+£t) 



\2S,i) 



(12) 



Eq is the energy of a pure 15 condensate with N — q atoms (see Eq. and |A5| ), Si = {i\ §^ + ^23'^ i^) K)' ^^"^ 
effective potential for the 25 atoms is 



y,^//(r).nr) + ^^^^^n.-,(r). 



(13) 



The density of 15 atoms remaining is n7v-g(r) = {N — q)\ip{r)\'' 
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Finding the 25* motional states which minimize with the requirement that they form an orthonormal 

basis, is equivalent to finding the eigenstates of the effective 25* Hamiltonian 



H. 



eff 
25 



(14) 



and the eigenvalue for state i is Si. The effective Hamiltonian (Eq. ^4|) is consistent with the two-component Hartree- 
Fock equations used to calculate the single particle wavefunctions for double condensates ||2l[] . The effective potential 
and some 25* motional states are depicted in Fig. 0. 

as Eq,i, using Eq. O and O, the energy supplied by two photons to 



If we denote the minimum of 
drive the transition to state i, for g <C A^, is 



2hv 



E. 



q,l 



En 



qiEis~2S + Si) + Eq - Eq 



q 

EiS~2S 



(15) 



We have used {Eq — -Eo)/? ~ dEo/dN = fi for small q. Note that < for states bound in the BEC interaction well. 
Since many 25* motional levels may be excited, there will be a distribution of excitation energies in the spectrum. 

When condensate atoms are coherently excited to an isolated level \i) by a laser pulse of duration t, the single 
particle wave functions evolve according to 



where 



Its', 0) ^ cos6i|lS', 0) +sin6'|2S',i) 



\{i\n{r)\0)\^ + Suj^ 



(|(^|^7(r)|0)| 



(16) 



(17) 



The detuning from resonance is Suj. In Eq. g^, we assume the excitation is weak enough to neglect the change in the 
single particle wave function for atoms in the condensate Depending upon which excitation scheme is being 

described, ri(r) is either f7Di?(r) or flDsir). 

The state vector for the system after excitation can be written 



\'i'{q),^) = (cos6'|lS',0) +sin6l|2S',i)) 



(cos6i|15, 0) +sin6l|2S',i)) 



N terms 



E 

9=0 



cos^-«0sin'?6 



m 



q\{N-q) 



(18) 



where the label (q) ~ Nsin^O is the expectation value of the number of 25 atoms excited. Although q is not a good 
quantum number for jvE'^g^ i), the spread in q, given by a binomial distribution, is strongly peaked around (q). 

For short excitation times, the population in state i grows coherently as t^. For the hydrogen experiment, however, 
although the excitation is weak and |(i|Jl(r)|0)| t ^ 1, i is longer than the coherence time of the laser (~ 200 /xs). 
This implies that the number of atoms excited to level i must be expressed in a form reminiscent of Fermi's Golden 
Rule. Equation n7\ can be rewritten in terms of a delta function using the relation sin^(xt)/7rx^t — > 6{x) as t — *■ cxd. 
(One can neglect^i|ri(r)|0) | compared to Suj because | (z|ri(r)|0) | is small compared to the spread in frequency of the 
laser excitation.) Then 



(q) ~ ^^mniv)\0)\^ 5{2hv ~ Eis^2S - + 



(19) 



It is understood that Eq. |l9| is to be convolved with the laser spectrum or a density of states function. The total 25* 
excitation rate is 



S{2hv) 



^ |(i|r!(r)|0)|2 5{2hi^ - Eis-2S - e^ + m) 
2 J2F'6i2hl^~Els-2S~£^+^J-)■ 
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(20) 
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Equation ^ defines the overlap factors, = |(i|f2(r)/r2 0)p, which are analogous to Franck-Condon factors in 
molecular spectroscopy. An expression equivalent to Eq. 2C|, the strength distribution function or dynamic form 
factor, is commonly used to describe collective excitations of many body systems Q. 

The BEC spectrum now appears as N times the spectrum of a single particle in |0) excited to eigenstates of the 
effective 25* potential. The broadening in the 15-25 BEC spectrum is homogeneous because it results from a spread 
in the energy of possible excited states, not from a spread in the energy of initially occupied states. 

The central results of this calculation arc the effective 25 potential (Eq. [l^) and the Fermi's Golden Rule expression 
for the excitation rate (Eq. po| ). Using this formalism we can now calculate the observed spectrum for Doppler-free 
and Doppler-sensitive excitation. 



D. Doppler-Free 1S-2S Spectrum 

Doppler-free excitation populates states which are bound inside the BEC potential well (see Fig. For a condensate 
in a harmonic trap, these states are approximately eigenstates of a three dimensional harmonic oscillator with trap 
frequencies larger than those of the magnetic trap alone by a factor of \/l — ais-2S /O'lS-is ~ 5 (see Eq. ^ and ^3|) . 
Because we know the wave functions, we can numerically evaluate Eq. gO[ The result of such a calculation is shown 
in Fig. |. 
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FIG. 2. Calculated Doppler-free spectrum of a condensate at T = in a three-dimensional harmonic trap. Zero de- 
tuning is the unperturbed Doppler-free transition frequency. The stick spectrum results from the sum over the transition 
amplitudes expressed in Eq. bd using the Thomas-Fermi density distribution for a peak condensate density of 10^^ cm~^ 
(47r?i^(ais-2S — ais-is)n{0)/m ~ 2/i x —0.95 MHz). The trap is spherically symmetric with ujtrap — 2tv x 6 kHz. The stick 
heights represent the coefficients of delta functions which must be convolved with the laser spectrum of about 1 kHz FWHM. 
The dashed curve (Eq. ^) follows from the integral over the BEC density distribution, Eq. |24| for the same peak condensate 
density. The envelope is independent of the symmetry of the trap, but the stick spectrum blends into a continuum in a trap 
with one weak confinement axis such as the MIT hydrogen trap [1,2]. Resolution of the individual transitions would require 
a stiff, near spherically symmetric trap, very stable experimental conditions, and high signal/noise. It does not seem feasible 
with the hydrogen experiment in the near future. 

At large red detuning {2h5v w A'KTT?ais-2sn{Q)/m — fi) transitions are to the lowest state in the BEC interaction 
well. The spectrum does not extend to the blue of 2h5v — because states outside the well have negligible overlap 
with the condensate and are inaccessible by laser excitation. In the overlap integrals in Fig. ^, wave functions for an 
infinite harmonic trap were used for the 25 motional states. These deviate from the actual motional states near the 
top of the BEC interaction well, introducing small errors in the stick spectrum nearer zero detuning. 

The envelope of the spectrum in Fig. ^ can be derived analytically and reveals some interesting physics. The 25 
single particle wave functions (r|i) oscillate rapidly. Thus the transition intensity to state i, governed by the overlap 
factor F}-,p = |(i|0)p, is most sensitive to the value of %l>{v) = (r|0) = y/n{r)/N at the state's classical turning points. 
At a given laser frequency, the excitation is resonant with all states with motional energy e = 2hv — Eis-2S- This 
suggests the excitation rate is proportional to the integral of the condensate density in a shell at the equipotential 
surface defined by the classical turning points of 25 states with motional energy e. 
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For a spherically symmetric trap, we can formally show this by making WKB and static phase approximations 
p5| , ^ - a technique which has recently been applied to describe s-wave collision photoassociation spectra |27| and 
quasiparticle excitation in a condensate One uses a WKB expression for the 25" eigenstate. Then, because of the 
slow spatial variation of the condensate wave function, the Doppler-free overlap factor only depends on the condensate 
wavefunction and the 15 and 25 potentials where the phase of the upper state is stationary. This yields 



DF 



47r 



Ri 



n{R,) 



N 



(21) 



where Ri is the Condon point, or the radius where the local wave vector of the excited state (fc2S = \l '^[^i — y2S^ (^)] ) 
vanishes. Ri is equivalent to the classical turning point for state i, and is defined through 



e, = V^^f{R.). 

Also, in the limit that we can neglect the slow spatial variation of the EEC wave function, D 

V^'s^iRi) is the slope of the effective 25 potential at the Condon point. 

Using Eq. |2l|, the Fermi's Golden Rule expression for the spectrum (Eq. pO) becomes 



(22) 



dV^^^{r)/dr\ 



SnFi^hiy) 



An 



Ri 



N 



■ 5{2hv - Eis-2S -£i + m)- 



(23) 



The Doppler-free excitation field and the BEC wave function are spherically symmetric, so only 25 motional states 
with zero angular momentum are excited. This implies that in the limit of closely spaced levels, — > J de 

in Eq. |2^. Using Eq. ^ we can change variables: J de = J dRV^g^^ (R) and 5 {2hv — Eis-2S — £ 

5 {2hv - Eis-2S - lZ!^Li£I]lM^ ^ where 5a = ais-2S - ais-is- This yields 



SoFi'ihv) = ""^^^DF j 47rdr r2n(r) ,5 {2hv - Eis-2S 



Anh^Sa n{r) 



(24) 



Using the probabilistic interpretation of |V'(ri)P (Sec. 11 B ), one can interpret Eq. gijin the following way. When a 
25 excitation is detected at a given frequency, it records the fact that a 15 atom was found at a position which had 
a 15 density which brought that atom into resonance with the laser. The rate of excitation is proportional to the 
probability of finding a condensate atom in a region with the correct density. This is a local density description of 
the spectrum, and it is justified by the slow spatial variation of the condensate wave function. 

For a Thomas-Fermi wave function in a three dimensional harmonic trap, Eq. e4 reduces to 



SoFi^hv) = 



IBirhnlpN {Eis- 



23 



(2h Sl'n 



2hv - Eis-2S 

2h5Vrn„,r. 



1/2 



(25) 



for 2h5vmax < '2'hv — i?is-2S < 0, and otherwise S£ip{2hv) = 0. Here, 2h5v„iax — Anh? Sa n{0) / m. 

Figure |^ shows that for a spherically symmetric trap, Eq. ^ agrees with the spectrum calculated directly with 
Fermi's Golden Rule (Eq. po[ ) using simple harmonic oscillator wave functions. For a trap which has a weak confinement 
axis, such as the MIT hyd rog en trap |jl|,||, discrete transitions in the spectrum are too closely spaced to be resolved. 
The envelope given by Eq. however, shows no dependence on the trap frequencies or the symmetry (or lack thereof) 
of the harmonic trap. 
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FIG. 3. Doppler-free spectrum of a condensate: comparison of theory and experiment (from [1]). The narrow feature near 
zero detuning is the spectral contribution from the noncondensed atoms (shown xl/40). The broad feature is the spectrum 
of the condensate. The dashed curve is Eq. |25| which comes from the integral over the BEG density distribution (Eq. p4| ) 
assuming a Thomas-Fermi density distribution for a harmonic trap. 

Theory and experimental data are compared in Fig. ^. Although the statistical error bars for the data are large 
due to the small number of counted photons, the theoretical BEC spectrum for a condensate at T = fits the 
data reasonably well. The deviations may indicate nonzero temperature effects or reflect experimental noise. The 
smoothing of the cutoff at large detuning may be due to shot to shot variation in the peak condensate density for the 
10 atom trapping cycles which contribute to this composite spectrum. Also, at low detuning the BEC spectrum is 
affected by the wing of the Doppler-frcc line for the noncondensed atoms. 

Using this theory, from the peak shift in the spectrum, the trap oscillation frequencies, and knowledge of ais-is 
and ais_2Sj one can calculate the number of atoms in the condensate. Assuming the experimental value of ais-2Sj 
the result is larger than the number determined from a model of the BEC lifetime and loss rates, which is discussed 
in [ p9[ . The uncertainties are large for these results, but the disagreement could be due to error in the experimental 
value of ais-2S: uncertainty in the gas temperature or trap and laser parameters, or thermodynamic conditions in 
the trapped gas which are different than assumed by the theories. For example, we have implicitly assumed local 
spatial coherence (g^'^HO) = I) [|o) in our form of the BEC wave function (Eq. [7|) . It has not yet been experimentally 
verified that the hydrogen condensate is coherent. 



E. Doppler-Sensitive 1S-2S Spectrum 



In contrast to the Doppler-free excitation spectrum, the Doppler-sensitive spectrum in principle reflects the finite 
momentum spread in the condensate as well as the mean field effects. The relevant momentum spread is given by 
the uncertainty principle and is ~ h/6z where 6z ^ 5 mm is the length of the condensate along the laser propagation 
axis. However, in the hydrogen experiment the cold collision frequency shift (~ I MHz) dominates over the Doppler- 
broadening in the spectrum {hko/2TnnSz « 100 Hz.) We can thus neglect Doppler-broadening, which is equivalent 
to neglecting the spatial variation of the BEC wave function in any transition matrix elements. In this regime it 



is possible to modify the derivation of the WKB and static phase approximations [g5-28| to calculate the Doppler- 
sensitive spectrum. 

We rewrite the Doppler-sensitive Rabi frequency (Eq. ^) as 

= 2nDs V4^(2?TT)z'ji(fcoO>l'"=°(^,0), (26) 

/ even 

where ji{kQr) is the spherical Bessel function of order I, and Yi"^{9,(p) is a spherical harmonic. This shows that the 
Doppler-sensitive laser Hamiltonian can excite atoms to 2S motional states with any even value of angular momentum, 
but with m = 0. 
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Transitions are to levels with motional energy ~ h^kQ/2m above the bottom of the 25* potential, so we label levels 
by A, their energy deviation from this value. For simplicity, we consider a spherically symmetric trap. This allows us 
to write a general expression for the 25* wave functions Vaj — Y"^=^{9, 4')u\^i{r)/r where UA,i{r)/r satisfies 



I) 



2m dr^ 2m 
Using Eq. |20|, the spectrum is 

NttTi 



UA.i{r) = Eis-2S + -TT^ + A UA.i{r) 
\ Zm J 



SDs{2hv) 



5Z l(^A,i|f^Ds(r)|0)p 5[2hu- Eis-2S 



A,l 



2m 



(27) 



(28) 



Using Eq. 



the overlap integral we must evaluate is 



(29) 



for I even, and otherwise. Because -y/ n{r) varies slowly, one can find an approximate expression for this matrix 
element. Appendix]^ gives the details of this derivation and uses the result to reformulate Eq. 28 as 

n(i?A) 



A ^2S l^Aj 



S{2hv - Eis-2S - 



2m 



A + ^i). 



(30) 



The matrix element (Eq. ^) gets it's main contribution at i?A where the classical wave vector of the WKB 
approximation for ma,z equals the classical wave vector of the WKB approximation for ji. In effect, i?A is the point 
where the spatial period of the wave function matches the wavelength of the laser field, 27r/fco (see Fig. This leads 
to a definition for i?A 



A = F;//(i?A), 



(31) 



which is identical to Eq. the definition of the Condon point from the calculation of the Doppler-free spectrum. 
Because the transition is localized in this way, the matrix element (Eq. |29| ) is proportional to \fn{RA)^ as is evident 
in Eq. M 
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FIG. 4. Effective potentials, wave functions and tfie laser field for Doppler-sensitive excitation of condensate atoms. Tlie 
spatial period of the 2S wave function, the laser wavelength, and the vertical axes for the potentials are not to scale. The 
vertical axes for the wave functions and laser field are arbitrary. In the overlap integral for the transition matrix element 
(Eq. ^9]), the only nonzero contribution comes from the region where the spatial period of the 25* wave function matches the 
wavelength of the laser field. This is indicated by the locations of the light vertical lines. As the laser frequency is changed, 
the region of wavelength match moves. 

Using Eq. we can replace the sum in Eq. |3^ with an integral and change variables, Sa J dA = J dR V^'g''^ {R). 
This yields the Doppler-sensitive lineshape 

Sosi^hi^) ~ TrUnls f ^^dr r^n{r) S ( 2h,, ~ E,s-2S - ^ - ^Man^\ ^^^^ 
J \ 2m m J 

The Doppler-sensitive condensate spectrum has the same shape as the Doppler-free spectrum, but it is shifted to the 
blue by photon momentum-recoil. Because ^ds — ^df/^, the Doppler-sensitive spectrum is half as intense as the 
Doppler-free. 

In pa] , experimental data are compared with Eq. p3, and the agreement is good. 



III. OTHER APPLICATIONS OF THE FORMALISM 
A. Other Atomic Systems and Excitation Schemes 

We have specifically considered 15-25 spectroscopy of hydrogen, but the formalism is more general. For instance, 
if the ground-excited state interaction were repulsive, this would simply modify the effective 25 potential (Eq. [T^ ) 
and the form of the motional states excited by the laser would change. Equations ^ and ^ would still be accurate 
for two-photon excitation to a different electronic state when the mean field interaction dominates the spectrum. 

In the recently observed rf hyperfine spectrum of a rubidium condensate pl| ] , the lineshape is determined by mean 
field energy and the different magnetic potentials felt by atoms in the initial and final states. The theory presented 
here can be modified to describe this situation as well. 

For Bragg diffraction or spectroscopy as performed in p^j3^ , atoms remain in the same internal state after exci- 
tation. Particle exchange symmetry of the wave function modifies the mean field interaction energy of the excited 
atoms with the atoms remaining in the condensate. In terms of the hydrogen levels, 15, F = 1, m/ = 1 atoms 
not in the condensate experience a potential of 8TTh'^ais-isnis(j^)/m. This is to be compared with the mean field 
potential of 47r?i^ai5_25?^is(r)/TO experienced by 25 particles excited out of the condensate and AT:h'^ais-isnis(j)/'m 
experienced by 15 atoms in the condensate. In appendix the point in the derivation where the difference arises is 
indicated. 



B. Doppler Broadening in the Doppler-Sensitive Spectrum 

To derive the Doppler-sensitive 15-25 spectrum, we neglected the variation of the condensate wave function, which 
is equivalent to neglecting the atomic momentum spread. This is well justified for the hydrogen experiment. The 
effect of small but nonnegligible momentum is discussed at the end of appendix |^. Now we briefly describe the 
Doppler-sensitive lineshape when Doppler- broadening is dominant. The lineshape turns out to be similar to that 
which was seen with Bragg spectroscopy of a Na condensate . 

When the mean field potential can be neglected, the 25 motional wave functions are approximately those of the 
simple harmonic oscillator potential produced by the magnetic trap alone. Because the spatial extent for these 
motional states is large compared to Sz^ in the region of the condensate the wave functions can be represented as 
plane waves momentum eigenstates |p) |34| . The spectrum becomes 

5z55(2H « ^ ^ |(p|r!z,5(r)|0)p S (2hi, - E,s-2S - |^ + m) 

^o-^ / '^'P l^(P - ^koi)\'S (2hv - E,s-2S - 1^ + 



p 



(33) 



The Fourier transform of the condensate wave function, yl(p) = / d^r e-''P-""/''7/'Ar(r), is nonzero for jp^;! <,?i/|(5a;|, 
\Pv\^^/\Sy\, and |p^| <7i/|(5z|. 
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The excited states have Pz « hko, so we define Sp = Pz — Tik^. Because the laser wavelength is smah compared to 
the spatial extent of the condensate, p^ /2m w ?i^fcg/2TO + hkoSp/m and the spectrum reduces to 

S{2hu) « ^^^^ J dp:.dpy\A{p,ic + Pyy + Sp{i.)z)\\ (34) 



where 



= 2hv - E,s-2S - ^ + M (35) 

TO ZTO 

defines the momentum class that is Doppler shifted into resonance. The spectrum is centered at 2hiy = -Eis-2S + 
h'^kQ/2m — fj,, and the lineshape depends on the orientation of the condensate wave function with respect to the laser 
propagation axis. 

For a Thomas-Fermi wave function in a spherically symmetric harmonic trap |A(p)p ~ b2(p''o/^)/(p^o/^)^P 



[ pO| , where tq — y2n{0)U/mw^. Numerical evaluation of the integral over p^ and Py shows that the lineshape is 
approximately given by the power spectrum of the wave function's spatial variation along z, S{2hv) cx |A((5p(i/)z)p. 

In recent experiments with small angle light scattering [^5|] , the momentum imparted to atoms is small compared 
to ^/2mCs^ where Cs — \J ii/m is the speed of Bogoliubov sound. In this case one can excite quasiparticles in the 
condensate as opposed to free particles. The theory described in this article only treats free particle excitation, but 
Bogoliubov formalism, combined with WKB and static phase approximations, has been used to describe the spectrum 
for quasiparticle excitation [ p8[ . 

IV. DISCUSSION 

To make the problem analytically tractable, we have only derived the BEC spectrum for the specific case of a 
spherically symmetric trap. The trap shape does not appear in the final expressions (Eq. ^ and |3^), however, and 
with reasonable confidence we can extend the results to any geometry. In the experiment, the trap aspect ratio is 
as large as 400 to 1, but the data agrees well with this theory. The physical picture of the transition occurring at 
the classical turning points, and the probabilistic or local density interpretation of the spectrum also support the 
generalization of Eq. ^ and ^ to 

5Di=^(2H = ^ / d3rnis(r)(5 \2hv - E^s~2S m ) ^ ^ 

Sus{2hv) = nhnls f d'rn,sir)S (2hu - E,s-2S - ^ - ^nh'San^s{r) \ 
J \ 2m TO / 

Equations |3^ and ^ take iTrfi^ 6a nis{r) /m as a local shift of the transition frequency and ascribe the excitation 
to a small region in space where the laser is resonant. This approach is similar to a quasistatic approximation 
in standard spectral lineshape theory [p6| which neglects the atomic motion and averages over the distribution of 
interparticle spacings to find the spectrum. Atom pairs at different separations experience different frequency shifts 
due to atom-atom interactions. This broadens the line. 

There are important differences between the theory presented here and the quasistatic approximation, however. 
For the standard quasistatic treatment to be valid, the lifetime of the excited state should be shorter than a collision 
time For a condensate, the classical concept of a collision time is inapplicable. We have shown that Eq. ^ and 
|37| result from a different approximation: neglecting the slow spatial variation of the BEC wave function. Also, for the 
condensate spectrum, one integrates over atom position in the effective potential, as opposed to integrating over the 
distribution of atom- atom separations. Finally, the BEC spectral broadening is homogeneous, which is not normally 
the case when making the quasistatic approximation. 

It is interesting that although the atoms in the condensate are delocalized over a region in which the density varies 
from it's maximum value to zero, the rapid oscillation of the excited state wave function essentially localizes the 



transition (Eq. 22 and p39| ). In this way, the excitation probes the condensate wave function spatially. 

The description of the BEC spectrum developed here has provided insight into the excitation process and it is 
general. We have shown that the formalism of transitions between bound states of the effective potentials can be used 
when either the mean field or Doppler broadening dominates. It can describe a variety of excitation schemes such as 
two-photon Doppler-free or Doppler sensitive spectroscopy to an excited electronic state, or Bragg diffraction which 
leaves the atom in the ground state. 
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APPENDIX A: ENERGY FUNCTIONAL FOR THE SYSTEM AFTER LASER EXCITATION 



In this appendix we derive Eq. the energy functional for the system after excitation which is minimized to find 
the 25* wave functions. 

The Hamiltonian and the excited state vector, are defined in Eq. |l| and |ll|. The symmetry operator is 

9 



explicitly written as S 



P, where the sum runs over the ^ ^ ^ 



q\{N~qy. 



distinct particle label 



permutations P. The energy functional for N — q IS condensate atoms and q 2S atoms in state i is 

N / 7 \ 



P 



(7V-g)(15,0| ( |- + ^(r) + i/-* 



1 15,0) 



^ '25,z) + ($,,,|i7^°"|$,,,). 



We evaluate the interaction term. 



(25;. 
(25;. 

(25,. 

(25;. 



15 
15 

15: 

15 



|5iJ^°"5|25;...;15;...) 
.|iJ^°"55|25;...;15;...) 

.|i/™"5|25;...;15;...)(^) 

.|iJ-"5]P|25;...;15;...), 



1/2 



(Al) 



(A2) 



where we have used [i/™",5] = and 55 ^ j''". Of the 7V(7V-l)/2 terms in i?'^"" (Eq. |), {N -q){N -q-l) /2 

of them result in a 15-15 interaction, {N — q)q of them result in a 15-25 interaction, and the rest result in a 25-25 
interaction which we can neglect. For the 15-15 terms, only the identity permutation contributes. For the 15-25 
terms two permutations contribute - the identity and switching the labels on the two interacting particles. The 
expectation value of H'^"''^ thus reduces to 



1/2 



{N - q){N - q - l)ais-is(0; 0\6{n - r2)|0; 0) + q{N - 9)ais-2s(»; 0|(5(ri - r2)|z; 0) 



(A3) 



As mentioned in Sec. Ill A. Eq. A2 would be modified for Bragg diffraction or spectroscopy as performed in p2 33 
because the internal state is unchanged during laser excitation. We do not explicitly treat this situation because it is 
not central to this study. 



Inserting Eq. A3 into Eq. Al, we find the energy functional is 
($,,.|il|<f,..) = E'o 



+q (2S,i\ 



p 



2m 

= E'a + q{EiS-2S + £i), 



V{v) 



UN-qir) 



\2S,t) 



(A4) 



where 
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E',^iN~ q){lS, 0| (f- + V{v) + W-') \1S, 0) 

+ ^{N ^q){N~q- l)ais-is(0; 0|5(ri - r2)|0; 0), 



(A5) 



is the energy ior N — q isolated 15 condensate atoms, and Ei 



2m ^2s'^ i^) I*)- "^^^^ density in the condensate 



for — g condensate atoms (g ^ N) is nN-q{r) — (N — g)(0|(5(ri — r)|0). 



APPENDIX B: WKB AND STATIC PHASE APPROXIMATIONS FOR THE DOPPLER-SENSITIVE 

EEC SPECTRUM 



In this appendix we calculate the Doppler-sensitive overlap integral, Eq. and simplify Eq. The derivation is 
similar to the treatment of p^j2^ . 

The overlap integral we must evaluate is 



Ia,i = (VA,i|e''"'+c"''"lO) 



(Bl) 



for / even and otherwise. 

Because ua,i and ji are rapidly varying compared to it is useful to express ma,; and ji in phase-amplitude form 
through a WKB approximation. We define the local wave vectors for ua.; and ji 



kj{ko,l,r) 

Then, in the classically allowed region 



P - 



^g + l) _ 2m ( ff 



"■0 ^2 
,2 



1/2 



1/2 



2m 



1/2 



sin/3„(A, Z,r), 



where 



/3„(A,/,r) = / dr'A:„(A,;,r')- V4, 



f3j{koJ,r) 



dr' kj{kQ, I, r') — 7r/4 



(B2) 
(B3) 

(B4) 
(B5) 

(B6) 
(B7) 



are the phases. The inner turning points against the centrifugal barriers are denoted by Rt- Note that the approxi- 
mations are good for (fcor)^ > l{l + 1). For (fegr)^ < l{l + 1), neglecting the small V^g'^ and A, the functions behave 
as damped exponentials. The outer turning points are of no concern to the calculation. 
Now we write 



lA,leven = -2^/4tt{21 + 1) / dr 



1 /2 

rt(r) sin/3i,( A, Z, r) / 2m \ sin/3j(/co, ^, J') 



V47r(2/ -f 1) ( ^ 



1/2 



dr , 



n(r) 



Nku{A,l,r)kokj{ko,l,r) 



cos [/3„(A,Z,r) -/3j(fco,;,r)]. 



(B8) 
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We have used the fact that \Jn{r) varies slowly and have dropped rapidly oscillating terms in the integral. 

We make the static phase approximation that the overlap integral will only have contributions from the point _Ra 
where the difference in the phase factors is stationary. This point is defined by = ^ (/?„ — I3j) — ku{A, I, Ra) — 
kj{ko, I, Ra), which is equivalent to an /-independent relation defining Ra for excitation to states with energy defect 
A, 



This is essentially identical to Eq. ^ from the calculation of the Doppler-free spectrum. 

We expand the difference in the phases in a Taylor series around Ra and write the overlap integral as 



(B9) 



V442^ + l)(^ 



1/2 



n(i?A) 



Nkok^ikoJ.RA) 



dx cos 



2h^kj{ko,l,RAy 



Pu{A,1,Ra)- PjikoJ,RA)- 



I 167r(2; + l)n(i?A) 
Nkok,{kQ,l,RA)V^s"iR^) 



-cos [/?„(A,Z,i?A) -/3i(fco,^,^A) ~7r/4] 



(BIO) 



To obtain the last line we have used the Fresnel integral dx cos(a + bx"^ 



^TT/bcos{a + TrT7r/4). Equation [BIC 



only holds for l{l + 1) < (fcoi?A)^. For l{l + 1) > (fcoi?A)^ Ia, 
Ra- 

From Eq. H and plO| , 



because ji{kor) is exponentially damped at 



A.l even 



167r(2; + l)n(i?A) 



2m 



2 A.l^en kok,{ko,l,RA)V^s"iRA) 

xcos2 [/3„(A, I, Ra) - /9, (fco, I, Ra) - tt/A] 6 (^2h,y - Eis~2S 
We can replace the cos^ function with it's average value of 1/2 because its phase varies rapidly with /. Thus 

(2^ + 1) 



(Bll) 



;(i+i)<(fcoflA) 

E 

I even 



kokj{ko,l,RA) 

/(i + l) = (fcoflA)^ 



cos^ [/3„(A,/,i?A)-/3j(fco,^,i?A)] 
1) 



klJl 



1(1+1) 



and 



i?A/2, 



2 1,2 



25 



n{RA) „ h^ka 

b[2hv - Eis-2S 



K7'{Ra) 



2m 



A + m). 



(B12) 



(B13) 



In the derivation given above, we neglected the variation of the condensate wave function, which is equivalent to 
neglecting the atomic momentum spread ~ h/5r, where 5r is the r extent of the condensate. When mean field effects 
dominate the spectrum, but the atomic momentum is not completely negligible, the lineshape wi ll deviate from Eq. 
^ only for small detunings, Si'^fikQ/2Trm Sr. One can see this from the overlap integral (Eq. B8) by expressing the 
condensate wave function in terms of the radial Fourier components, Ar{p) — J dr e~^P^^^'ip(r)^ to obtain 



Ia 



I even 



-2^47r(2/ + l) 



2'Kh 



dpMp)ldre^^^'^'^^^dtM(^ 



sin/3j-(fco,/,r) 
y/kokj{ko,l,r) 



(B14) 
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Each momentum compone nt wi ll only contribute to the matrix element at the point Ra,i,p where the total phase 
under the r integral in Eq. |B14 is stationary. This leads to a definition of Ra,i,p for each momentum, p/h = 
\ku{A,l, Ra,i,p) — kj{kQ,l, RAd,p)\- When |A| ^ h^ko/mSr, p/h is negligible and this yields the same relation as 
found by neglecting the curvature of the BEC wave function (Eq. p9| ). This implies S{2h\6i'\ ^ fi^ka/mSr) is 
unaffected by the atomic momentum. When \l!^\<Ti'^kQ/m5r^ the momentum spread in the condensate alters Ia,i even- 
Thus S{2h\5v\<h^ko/m5r) will show some Doppler-broadening because of finite atomic momentum. This effect is 
negligible for the hydrogen condensate because the cold collision frequency shift 1 MHz) is much greater than the 
Dopplcr width resulting from a 5 mm long condensate wave function (?ifco/27rm(5z ~ 100 Hz). 
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